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Abstract 

Bose-Einstein correlations in one and two dimensions have been studied in deep 
inelastic ep scattering events measured with the ZEUS detector at HERA using 
an integrated luminosity of 121 pb^^. The correlations are independent of the 
virtuality of the exchanged photon, Q^, in the range 0.1 < < 8000 GeV^. 
There is no significant difference between the correlations in the current and 
target regions of the Breit frame for > 100 GeV^. The two-dimensional shape 
of the particle-production source was investigated, and a significant difference 
between the transverse and the longitudinal dimensions of the source is observed. 
This difference also shows no dependence. The results demonstrate that 
Bose-Einstein interference, and hence the size of the particle-production source, 
is insensitive to the hard subprocess. 
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1 Introduction 



The quantum-mechanical wave-function for a pair of identical bosons is symmetric under 
particle exchange. As a consequence, interference effects are expected between identical 
bosons emitted close to one another in phase space. These effects enhance the two-particle 
density at small phase-space separations. Such Bose-Einstein (BE) correlations for like- 
charged hadrons were first observed by Goldhaber et al. P in pp annihilation. 

The BE correlations in momentum space are related to the spatial dimensions of the 
production source. Therefore, studies of the BE effect may lead to a better understanding 
of the structure of the source of the identical bosons. In deep inelastic scattering (DIS), the 
production volume may depend on the virtuality of the exchanged photon, Q^, since the 
transverse size of the virtual photon decreases with increasing Q^. Under this hypothesis, 
the BE correlations would depend on Q^. 

Alternatively, the size of the region over which BE correlations take place may be de- 
termined only by the soft, or fragmentation, stage of the process. The BE effect is 
independent of if hard scattering and fragmentation factorise. For example, in the 
Lund fragmentation model [210]; no sensitivity to is expected, and the BE correla- 
tions between identical bosons are a measure of the tension of the colour string between 
partons. 

This paper reports on investigations of the BE correlations in one and two dimensions 
in neutral current ep DIS, focusing on the dependence of the BE interference on Q^. 
The BE effect in one dimension is measured with a higher precision than previously 
in ep collisions |3] and over a much wider kinematic range, from ~ 0.1 GeV^ to 
8000 GeV^. The correlations are also studied for the first time in DIS in the current and 
the target fragmentation regions of the Breit frame 0, which are known to have rather 
different fragmentation properties jH]. The two-dimensional analysis provides sensitivity 
to a possible elongation of the source expected in the Lund string model [3]. 

2 Definition of measured quantities and 
model predictions 

Bose-Einstein correlations are usually parameterised using a Gaussian expression for the 
normalised two-particle density \7\ : 

RiQ,2) =ail+P Qu) (1 + A e-^' (1) 

where Q12 = a/— (pi — ^2)^ = a/ — 4:m'l^^^^ is the Lorentz-invariant momentum differ- 
ence between the two identical bosons, which is related to the invariant mass, M , of the 
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two particles with four-momenta pi and p2 and mass mboson- In the present analysis, all 
charged particles are assumed to be pions. The parameter A is a measure of the degree 
of coherence, i.e. the fraction of pairs of identical particles that undergo interference. 
For a totally coherent emission of pions, A = 0, while for an incoherent source, the sym- 
metrisation of the wave function for identical particles leads to A = 1. The quantity r is 
the radius of the production volume, while the phenomenological parameter /5 is used to 
take into account any long-distance correlations, and a is a normalisation constant. The 
Gaussian parameterisation in Eq. is motivated by the assumption that the emitting 
sources of identical bosons are described by a spherical Gaussian density function. 

When the BE correlations are interpreted in terms of the Lund fragmentation model, the 
correlation strength is related to the string tension. In this case, the correlations should 
have an approximately exponential shape i.e. \e~^^ '^12 in Eq. (0) should be replaced 
by X'e-''^'\ 

To calculate R^Qu), the inclusive two-particle density, p = 1/N ■ dnpairs/dQi2, was used, 
where npairs is the number of particle pairs and is the number of events. The densities 

were calculated for like-charged particle combinations {p{++, )) and for unlike-charged 

combinations (p(H — )), and the ratio computed as = p(++, )/p(H — )• This technique 

helps to remove correlations due to the topology and global properties of the events 

contributing to ). The quantity C, contains additional short-range correlations, 

due to resonance decays (contributing to p{-\ — )), which must also be removed. To reduce 
such non-BE effects, a Monte Carlo (MC) sample without the BE effect was used to 
calculate ^'^c.noBE^ ^^^^ non-Bose-Einstein effects were removed by use of the double ratio, 
R{Qi2) = ^<i^t^/^^^c,noBE_ ^]^jg casc, it was assumed that the BE effect can be factorised 
from other types of correlations, and that non-BE correlations are well described by the 
MC models. 

It is also possible to extract the BE parameters by considering pairs of particles from 
different events as a reference sample (track-mixing method). This method was not used 
in the current analysis since it is difficult to control the systematic effects arising from 
differences of the events in the sample. Earlier experiments jHMl] have found that values of 
r extracted with the track-mixing method is systematically smaller than those obtained 
with the method used here. 

Bose-Einstein correlations can be studied in two dimensions. For this, the longitudinally 
co-moving system (LCMS) jH] is often used, since in this frame the correlations have a 
convenient interpretation in the Lund-string model. For ep collisions, the LCMS can be 
defined for each pair of particles with momenta pi and p2 as the frame in which the 
sum of the two momenta, pi + ps, is perpendicular to the 7*g axis, as shown in Fig. ^ 
The three- momentum difference, Q = (p2 — Pi), can be decomposed in the LCMS into 
transverse, Qt, and longitudinal, Ql, components. The longitudinal direction is aligned 
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with the direction of motion of the initial parton. Therefore, in the string model, the 
LCMS is the local rest frame of a string. In this case, the BE effect can be parameterised 
using the two-dimensional function 

R{Qt, Ql) = « (1 + a Qt + a Ql) (1 + a e~'T QT~rl Qi)^ (2) 

where tt and are the transverse and longitudinal size of the boson source. The mea- 
surements were done using this parameterisation and the same procedure as for the one- 
dimensional study. 

The Lund model predicts that the longitudinal size of the production source is larger than 
the transverse one, t^, > r^- However, the usual implementation of the BE effect in 
the Lund MC models does not contain such an elongation. 

The Breit frame jH], as shown in Fig. Q allows the separation of the radiation of the 
outgoing struck quark from the proton remnant. Therefore, this frame can be used to 
test the sensitivity of the BE effect to different underlying dynamics. All particles with 
negative p^'''^'* form the current region. These particles are produced by the fragmentation 
of the struck quark, so that this region is analogous to a single hemisphere of an e^e~ 
annihilation event, while the target region is dominated by the softer fragmentation of 
the proton remnants. 



3 Experimental setup 

A detailed description of the ZEUS detector can be found elsewhere ^U]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) ^T], which operates 
in a magnetic field of L43T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift chamber layers, organized in nine superlayers covering the 
polar- angle^ region 15° < ^ < 164°. The transverse-momentum resolution for full-length 
tracks is (t{pt)/pt = 0.0058pT © 0.0065 © 0.0014/pt, with pr in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) J2] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a{E)/E = 0.18/ \/E for electrons 
and a{E)/E = 0.35//E for hadrons, with E in GeV. 

A presampler ^2] is mounted in front of FCAL, BCAL and RCAL. It consists of scin- 
tillator tiles that detect particles originating from showers in the material between the 
interaction point and the calorimeter. This information was used to correct the energy 
of the scattered electron. The position of electrons scattered close to the electron beam 
direction is determined by a scintillator strip detector (SRTD) |14j . 

The beam-pipe calorimeter (BPC) ^3] was installed 294 cm from the interaction point in 
order to enhance the acceptance of the ZEUS detector for low-Q^ events. It is a tungsten- 
scintillator sampling calorimeter with the front face located at Z = —293.7 cm, the centre 
at Y = 0.0, and the inner edge of the active area at X = 4.4 cm, as close as possible 
to the rear beam pipe. The energy resolution as determined in test-beam measurements 
with 1-6 GeV electrons is as/E = 17%//E, with E in GeV. 

4 Data sample 

Two data samples were used for the present analysis. The data for > A GeV^ were 
taken during the 1996-2000 period and correspond to an integrated luminosity of 121 
pb~^. The lepton beam energy was 27.6 GeV and the proton beam energy was 820 GeV 
(1996-1997) or 920 GeV (1998-2000). The second sample consists of low-Q^ events taken 
with the BPC. This sample corresponds to 3.9 pb~^ taken during 1997. 

The kinematic variables and Bjorken x were reconstructed using the electron method 
(denoted by the subscript e), which uses measurements of the energy and angle of the 
scattered lepton, the double angle (DA) method [12] or the Jacquet-Blondel (JB) method 

[HI- 

The scattered-lepton candidate for the region > 4 GeV^ was identified from the pattern 
of energy deposits in the CAL ^Hj- The following requirements were imposed: 

• g2 > 4 GeV2; 

• E^r > 10 GeV, where E^' is the corrected energy of the scattered lepton measured in 
the CAL; 

• 40 < (5 < 60 GeV, where ^ = ~ cos^j), Ei is the energy of the ith calorimeter 
cell, 6i is its polar angle and the sum runs over all cells; 

• ye< 0.95 and yjB > 0.04; 

• I 2'vertcx I < 50 cm, where Z^^rtcx is the vertex position determined from the tracks; 
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• the position of the scattered-lepton candidate in the RCAL was required to be outside 
a box of ±14 cm in X and Y . 

In total, 6.4 milhon events were selected. 

The low-Q^ events, selected in the region 0.1 < Ql < 1.0 GeV^, were reconstructed by 
identifying energy deposits in BPC consistent with a scattered positron with an energy 
of least 7 GeV. The positron position at the BPC front face had to lie within the fiducial 
area, 5.2 < X < 9.3 cm and —2.3 < Y < 2.8 cm. Other cuts are identical to those used for 
the data sample at > 4 GeV^, except for the yjB cut, which was raised to i/jb > 0.06. 
The reconstruction of the Breit frame and the variable were performed using variables 
calculated with the electron method. In total, about 100000 events were selected. 

The measurement uses CTD tracks assigned to the primary event vertex. Tracks were 
required to pass through at least three CTD superlayers and have transverse momentum 
Py*^ > 150 MeV. The approximate pseudorapidity range for selected tracks is \7]\ < 1.75. 
To ensure that tracks were well reconstructed, track pairs were required to satisfy Q12 > 
0.05 GeV. In the kinematic region used for this analysis, pions constitute about 81% of 
the tracks, according to MC expectations. 

5 Reconstruction procedure and 
acceptance correction 

The parameters for the BE effect were determined from a fit using either the Gaussian or 
exponential parameterisation. The reference sample was calculated using unlike-charged 
pairs, and then the MC was used to remove the effect of resonances, as discussed in Sect. |21 
The regions affected by imperfections in the MC simulations of and p° decays were 
excluded from the fit. 

The measured correlation functions, R{Qi2) and R{Qt,Ql), were corrected for detector 
effects using a bin-by-bin procedure. The detector correction factor was calculated from 
MC events as ^sc^^^dct^ where ^s<='i(^'i<=*) is the generated (reconstructed) two-particle 

density ) divided by p(H — ). No cuts were applied on the true MC hadrons. 

The corrections are close to unity (1.0 ± 0.1), since some detector effects cancel in the 
ratios of the two-particle densities. Only for the lowest Q12 region (0.05 < Qu < 0.1 
GeV) the correction is as large 1.3-1.4. 

Since the generator-level MC events do not have the BE effect, and are used both for the 
acceptance correction and for the subtraction of resonances, the final result is equivalent 
to the detector-level measurement in the restricted kinematic region defined by the cuts 
in Sect. |3] The effect of these cuts, in particular those on track momentum and angle. 
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depends upon Q^, since determines the available phase space for particle production. 
The MC model containing BE correlations was used to estimate the effect of extrapolation 
to the full phase space available in each bin. The effect was found to be small (+3.5% 
for the BE radii) and independent. This model-dependent correction was not included 
in the final results. 

The MC events were generated with the ARIADNE 4.08 model ^H] interfaced with HER- 
ACLES 4.5.2 [201 using the DJANGOH 1.1 program [21] to incorporate first-order elec- 
troweak corrections. The generated events were then passed through a full simulation of 
the detector using GEANT 3.13 |22] and processed with the same reconstruction program 
as the data. The detector-level MC samples were then selected in the same way as the 
data. 

To simulate hadronisation, ARIADNE employs the PYTHIA 5.7/JETSET 7.4 program 
|23j . which is based on the Lund string model |23- The BE effect, which is treated as 
a final-state interaction by redistributing hadron momenta according to a chosen param- 
eterisation, is available as an option in the PYTHIA/ JETSET program. The default 
parameter setting for ARIADNE does not contain the BE effect. For systematic checks, 
the BE correlations were included in the simulation for the acceptance calculations. The 
BE effect was parameterised by a Gaussian function with the parameters as determined 
by HI U. 

As a systematic check, the HERWIG 6.2 |23] model was used to calculate the acceptance 
correction. The hadronisation stage in HERWIG is described by a cluster fragmentation 
model [2n]- The BE effect is not implemented in the HERWIG model. 

6 Systematic uncertainties 

The systematic uncertainties of the measured BE correlations were determined by chang- 
ing the selection cuts or the analysis procedure and repeating the extraction of the BE 
parameters. The following systematic studies have been carried out, with the resulting 
uncertainty for the BE radius for the highest-precision measurement (at > 4 GeV^) 
given in parentheses: 

• the event-selection cuts on ?/e, yjB-, -^vertex, Ee' were varied, and the scattered- lepton 
energy scale was changed within its uncertainty ±2% (^q4%). 

• the DA method was used to reconstruct and the boost vector to the Breit frame 
(-0.4%); 

• the minimum transverse momentum for tracks was raised by 50 MeV (—3%); 
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• tracks were required to be within the pseudorapidity range | |< 1.5, in addition to 
the requirement of three CTD superlayers (+0.1%); 

• the fit was performed excluding one data point from each side of the region affected by 
p decays. In addition, the lowest Q12 bin, which is not well described by the Gaussian 
parameterisation was excluded from the fit (1q 4%). 

The overall systematic uncertainty was determined by adding the above uncertainties in 
quadrature. 

As an additional systematic check, the HERWIG model was used to calculate the accep- 
tance correction. Change in the reconstructed BE radius was +9%. However, HERWIG 
does not reproduce the measured p(H — ) density, so the extraction of the BE effect is less 
reliable; therefore, this check was not included in the final systematic uncertainties. 

Identically charged particles are subject to Coulomb repulsion, which is not simulated by 
MC models. The Bose-Einstein correlation function was corrected in the data using the 
Gamow factor [2ZI. After correcting for the Coulomb effect, the size of the BE radius and 
A slightly increased (+4% for r). It was found that this correction does not depend on 
Q^, therefore, it was not included in the final results. 

7 Results and discussion 
7.1 One-dimensional study 

Figure 121 shows the measured R{Qi2) for Q"^ > 4 GeV^ together with the Gaussian fit of 
Eq. ((H) and the exponential fit. The observed distortions for Q12 > 0.9 GeV are caused 
by the decay products of resonances which are not well described by the MC simulation. 
Both parameterisations give fits of similar quality. The extracted parameters for the 
Gaussian and exponential fits are given in Table H 

The BE parameters extracted using the Gaussian parameterisation are shown in Fig. El as 
functions of Q^. Within the statistical and systematic uncertainties, the data indicate no 
variations with the virtuality of the exchanged photon in the range 0.1 < < 8000 GeV^ 
when either the Gaussian or exponential parameterisation is used. This measurement is 
consistent with an earlier HI measurement |1| for 6 < < 100 GeV^ using the wrong- 
charge background subtraction. 

Figure El also shows the comparisons between the current and the target regions of the 
Breit frame. The BE effect in the current region was extracted for > 100 GeV^, where 
the charged multiplicity is high enough for a reliable measurement. The result indicates 
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that there is no significant difference between the BE effects in the current and the target 
regions of the Breit frame. Note that the data shown in Fig. El for the total phase space 
are dominated by the target region. 

7.2 Two-dimensional study 

The BE correlations can be studied in more than one dimension after decomposing the 
momentum difference into its transverse and longitudinal components. The BE parame- 
ters, ri, tl and A are shown in Figure EJ Table |21 also gives the ratio tt/tl. The result 
shows that the pion-emitting region, as observed in the LCMS, is elongated with being 
larger than tt- 

Figure |3] and Table El also show extracted BE parameters as a function of for the 
two-dimensional measurements. There is again no dependence of the BE radii. 

7.3 Comparisons with other experiments 

One-dimensional BE correlations have been measured in a number of experiments using 
the wrong-charge background subtraction |281l29ll3UII311l321l33j . For DIS, the present result 
agrees well with the /ip data at > 4GeV^ measured by the EMC Collaboration |28j . 
as well as with the data for Q < < 100 GeV^ measured by HI |3], as discussed before. 

The data also agree with the LEPl average r = 0.78 ± 0.01(stat.) ±0.16{syst.) calcu- 
lated by combining results for the track-mixing and wrong-charge background-subtraction 
methods, as well as with LEP2 measurements [nDj. The present results also agree with 
earlier lower-energy e^e~ annihilation experiments (see a review jHSI)- Since the BE effect 
measured in this paper is dominated by the target region, this suggests that the BE effect 
is not sensitive to the underlying hard processes. 

Comparisons with 7T~^p and pp data [3211^ indicate that the BE radii may be larger for 
these two processes. However, since the observed BE interference significantly depends on 
the experimental procedure used to extract the effect, it is difficult to assess quantitative 
differences between the BE correlations observed in DIS and hadron-hadron collisions. 

The measured BE effect reported here disagrees with that in relativistic heavy-ion colli- 
sions, which are characterised by significantly larger BE radii, which depend on the atomic 
number. A, of the projectile as r ^ 0.7^1/3 fm m. 

The LEP experiments have recently reported an elongation of the pion source in e"'"e~ an- 
nihilation events [HBl- The present ratio for the two-dimensional BE effect in ep collisions 
is consistent with these measurements. 
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8 Conclusions 



One- and two-dimensional Bose-Einstein correlations have been studied in deep inelastic 
ep scattering. The effect was measured as a function of the photon virtuality, Q^, in the 
range from 0.1 to 8000 GeV^. The results indicate that the source of identical particles 
has an elongated shape, consistent with the expectations of the Lund model. 

The Bose-Einstein effect in one and two dimensions does not depend on the virtuality 
of the exchanged photon. The elongation of the pion source is also independent of Q^. 
In addition, the Bose-Einstein correlations in the current and target regions of the Breit 
frame are similar, even though there is a significant difference in the underlying physics 
in these two regions. 

These high-precision results, obtained over a wide kinematic range in a single experiment, 
demonstrate that Bose-Einstein interference in ep collisions, and hence the size of the 
particle-production source, is insensitive to the hard subprocess. 
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(GeV^) 
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r (rmj 


A 


r (imj 


A r\ r\r\ r\ 

4 - 8000 


0.475 ± 0.007„o;go3 


r\ 1 r\ r\r\r\~\~\\ \\/.}. 

0.666 ± 0.009_o:o36 


0.913 ± 0.015^°:°^^ 


0.928 ± 0.023lo;o94 


100 - 8000 


0.431 ± 0.015tJ]:[!J^ 


0.646 ± 0.02lt°:°°^ 


0.815 ± 0.037i[!:Ji^ 


0.859 ± 0.059t[!:?}^ 


0.1 - 1 


0.464 ± 0.027lJ5:[!^^ 


0.602 ± 0.036^!}:[!^? 


0.929 ± 0.069^^32 


0.785 ± 0.071+[!:J^^ 
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0.468 ± 0.020t[}:|][|^ 


0.685 ± 0.028^!}:[!^} 


0.892 ± 0.043^15;^^^ 


0.954 ± 0.069^^68 


8 - 16 
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0.620 ± o.oi8^!}:[!^^ 


0.911 ± 0.041+°:;^^! 


0.857 ± 0.0541°:°^° 


16 - 32 


0.473 ± 0.017^°:°J^ 


0.629 ± 0.022+°;°°^ 


0.926 ± 0.052+°:J[^ 


0.829 ± 0.066+°:?}^ 


32 - 64 


0.496 ± 0.018tJ]:|]?|j 


0.679 ± 0.022i°:°22 


0.941 ± 0.042+°:J°° 


0.910 ± o.o6ot[!:|]^^ 
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0.445 ± 0.0171^03 


0.665 ± 0.024+°:°°^ 


0.843 ± 0.038t[!:Jj^? 
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0.431 ± 0.021+°:°}° 


0.649 ± 0.030+!}:[!!|^ 


0.821 ± 0.050t[!:J^?^ 
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400 - 1200 
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0.889 ± 0.216+[!:J^^^ 
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Q2 (GeV^) 


A 


tl (fm) 


tt (fm) 




4 - 8000 


0.44 ± o.oilE5:[!3 


0.95 ± o.o3lE5:[!^ 


0.69 ± OMtofe 


0.72 ± 0.03lJ^:[!3 


100 - 8000 


0.32 ± 0.03+°:°^ 


0.88 ± 0.08+°:°^ 


0.62 ± 0.04+°:°? 


0.70 ± 0.08+°:°? 


0.1 - 1 


0.41 ± o.o5i[5:[!^ 


0.82 ± Omtom 


0.74 ± 0.08t[!:?^ 


0.91 ± 0.141°:°^ 


4 - 16 


0.46 ± 0.02+°:°? 


0.84 ± 0.04+°:°^ 


0.69 ± 0.02+°:°^ 
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16 - 64 


0.39 ± 0.02+°:°^ 
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400 - 8000 


0.42 ± 0.10+°:°? 


1.08 ± 0.27+°:^^ 


0.67 ± 0.1ll°:JJ 


0.62 ± 0.18+°:°^ 



Table 2: The values of the two-dimensional BE parameters for different 
ranges, as well as the ratio tt/tl. The Gaussian parameterisation was used to 
extract the parameters. The statistical and systematic uncertainties are indicated. 
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Figure 1: The longitudinally co-moving system for a pair of particles in DIS. 
This system is defined as the frame of reference in which the sum of the two-particle 
momenta, pi and p2, is perpendicular to the 'j*q axis, which coincides with the Z 
axis of the Breit frame. 
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Figure 2: The measured Bose-Einstein correlation function, R{Qi2), together 
with the Gaussian and the exponential fits. The error bars show the statistical 
uncertainties. The data points included in the fit are marked with the circles. 
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Figure 3: The extracted radius, r, and the incoherence parameter, X, as functions 
of for the total measured phase space (left figures) and for the target region and 
the current region of the Breit frame (right figures). The HI data point is shown 
for the mean value of the measured range, 6 < < 100 GeV^ . The BE effect 
is shown for the current fragmentation region for 100 < g2 < 8000 GeV^. The 
inner error bars are statistical uncertainties; the outer are statistical and systematic 
uncertainties added in quadrature. The dotted lines show the average values, with 
the x^/ndf ~ 0.5 for both r and X. 
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Figure 4: The extracted radii, tt, tl, and the incoherence parameter X as 
functions of for the two-dimensional correlation function R{Qt, Ql) ■ The in- 
ner error bars are statistical uncertainties; the outer are statistical and systematic 
uncertainties added in quadrature. The dotted lines show the average values. 
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